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Abstract 

The main challenge in the Standard Model calculation of the mass and width difference in the 
D° — D° system is to estimate the size of SU(3) breaking effects. We prove that D meson mixing 
occurs in the Standard Model only at second order in SU (3) violation. We consider the possibility 
that phase space effects may be the dominant source of SU(3) breaking. We find that y = Ar/2r 
of the order of one percent is natural in the Standard Model, potentially reducing the sensitivity 
to new physics of measurements of D meson mixing. 
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I. INTRODUCTION 



It is a common assertion that the Standard Model prediction for mixing in the D° — D° 
system is very small, making this process a sensitive probe of new physics. Two physical 
parameters that characterize D° — D° mixing are 

AM AT 

* = — , y = ^, (i) 

where AM and Ar are the mass and width differences of the two neutral D meson mass 
eigenstates, and V is their average width. The D° — D° system is unique among the neutral 
mesons in that it is the only one whose mixing proceeds via intermediate states with down- 
type quarks. The mixing is very slow in the Standard Model, because the third generation 
plays a negligible role due to the smallness of IV^V^ and the relative smallness of nib, and 
so the GIM cancellation is very effective |]TJ, 0, [| (|, 

The current experimental upper bounds on x and y are on the order of a few times 
10~ 2 , and are expected to improve significantly in the coming years. To regard a future 
discovery of nonzero x or y as a signal for new physics, we would need high confidence that 
the Standard Model predictions lie significantly below the present limits. As we will show, 
in the Standard Model x and y are generated only at second order in SU(3) breaking, so 
schematically 

x , y ~ sin 2 9 C x [SU(3) breaking] 2 , (2) 

where 8q is the Cabibbo angle. Therefore, predicting the Standard Model values of x and 
y depends crucially on estimating the size of SU(3) breaking. Although y is expected to 
be determined by Standard Model processes, its value nevertheless affects significantly the 
sensitivity to new physics of experimental analyses of D mixing ||. 

At present, there are three types of experiments which measure x and y. Each is actually 
sensitive to a combination of x and y, rather than to either quantity directly. First, there 
is the D° lifetime difference to CP even and CP odd final states |7], ||, §, [TU], |TTJ, which to 
leading order measures 

r(D -> n+K~) . A m 

y° p= T (D^K+K-) ~ I = ycos0-xsm0 — , (3) 

where A m = \q/p\ 2 — 1 (see Eq. ([]) for the definition of the neutral D mass eigenstates), and 
is a possible CP violating phase of the mixing amplitude. Second, one can measure the 
time dependence of doubly Cabibbo suppressed decays, such as D° — > K + n~ fl2]| , which is 
sensitive to the three quantities 

(x cos 5 + y sin S) cos (p , (y cos 5 — x sin S) sin <p , x 2 + y 2 , (4) 

where 5 is the strong phase between the Cabibbo allowed and doubly Cabibbo suppressed 
amplitudes. A similar study for D° — * K~tt + 7t° also would be valuable, with the strong 
phase difference extracted simultaneously from the Dalitz plot analysis |]13 |. Third, one can 



search for D mixing in semileptonic decays |14| , which is sensitive to x 2 + y 2 

In a large class of models, the best hope to discover new physics in D mixing is to observe 
the CP violating phase, 0i2 = arg[Mi 2 /ri 2 ] (see the definitions (|7|) and (||) below), which 
is very small in the Standard Model. However, if y 3> x, then the sensitivity of any physical 
observable to (fiu is suppressed, since A m is proportional to x/y and is to (x/y) 2 , even 
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if new physics makes a large contribution to AM |J. It is also clear from Eq. (Q) that if 
y is significantly larger than x, then S must be known very precisely for experiments to be 
sensitive to new physics in the terms linear in x and y. It may be possible to measure S 
with some accuracy at the planned r-charm factory CLEO-c [|1^, [T6|j . 

There is a vast literature on estimating x and y within and beyond the Standard Model; 
for a compilation of results, see Ref. [O. Roughly there are two approaches, neither of 



which give very reliable results because m c is in some sense intermediate between heavy and 
light. The "inclusive" approach is based on the operator product expansion (OPE). In the 
m c ^> A limit, where A is a scale characteristic of the strong interactions, AM and Ar can 
be expanded in terms of matrix elements of local operators fl], @; 5H- Such calculations yield 
x,y 10~ 3 . The use of the OPE relies on local quark-hadron duality, and on A/m c being 
small enough to allow a truncation of the series after the first few terms. The charm mass 
may not be large enough for these to be good approximations, especially for nonleptonic D 
decays. An observation of y of order 10~ 2 could be ascribed to a breakdown of the OPE 
or of duality ||I8 |, but such a large value of y is certainly not a generic prediction of OPE 



analyses. The "exclusive" approach sums over intermediate hadronic states, which may be 
modeled or fit to experimental data ||, [H| |20| . Since there are cancellations between states 



within a given SU(3) multiplet, one needs to know the contribution of each state with high 
precision. However, the D is not light enough that its decays are dominated by a few final 
states. In the absence of sufficiently precise data on many decay rates and on strong phases, 
one is forced to use some assumptions. While most studies find x, y < 10~ 3 , Refs. |22], [23| 



obtain x and y at the 10~ 2 level by arguing that SU (3) violation is actually of order unity, 
but the source of the large SU(3) breaking is not made explicit. 

In this paper, we compute the contribution to Ar from SU(3) breaking from final state 
phase space differences. This is a calculable source of SU(3) violation, which enhances the 
rates to final states containing fewer strange quarks. In Sec. we review the formalism of 
D° — D° mixing. In Sec. [TIT] we give a general group theory proof that AM and Ar are 
only generated at second order in SU (3) breaking if SU(3) violation enters these quantities 
perturbatively. In Sec. |TV| we discuss the estimates of SU (3) breaking using the "inclusive" 
and "exclusive" analyses, and remind the reader of the shortcomings of each. Our main 
results are found in Sec. [V|, namely the calculation of SU(3) breaking in Ar from phase 
space effects in two-, three- and four-body final states. We find that such effects are very 
important, and can naturally account for Ar/2r at the percent level. We extend the analysis 
to intermediate resonances in Sec. [V]]. In Sec. [VII| we present our conclusions and ask whether 
in light of our results it remains possible for the measurement of D mixing to probe new 
physics. 



II. FORMALISM 

We begin by reviewing the formalism for D° — D° mixing. The mass eigenstates and 
Ds are superpositions of the flavor eigenstates D° and D°, 

\D L , s )=p\D°)±q\W), (5) 

where \p\ 2 + \q\ 2 = 1. In the Standard Model CP violation in D mixing is negligible, as is 
CP violation in D decays both in the Standard Model and in most scenarios of new physics. 
From here on we will assume that CP is a good symmetry. Then p = q, and I-Dl.s) become 
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CP eigenstates, 



CP\D ± ) = ±\D ± ) 



(6) 



with the mass and width differences defined as AM = mo + — ttid_ and Ar = Td + — 
The off-diagonal element of the D° — D° mass matrix can be expressed as 



M 12 = (D \Ht C=2 \D°}+PY: 



(D°\Hu \n){n\H£; \D°) 



m 2 D - El 



T 12 = Y.Pnm< C=l \n){n\>Ht C = l \D») 



(7) 
(8) 



where the sum is over all intermediate states, P denotes the principal value, and p n is the 
density of the state n. The first term in Eq. (|?p comes from the local |AC| = 2 operators 
(box and dipenguin), which affect M 12 only. The second term comes from the insertion of 
two | AC | = 1 operators. There is a contribution of this type to both M\i and r^. 

One can then express y in two equivalent ways, either as a sum over the states that are 
common to D° and D°, 



y = ^Y,Pn [(D°\n w \n)(n\n w \D°) + (D°\n w \n)(n\n w \D°) 



(9) 



or as the difference in the decay rates of the two mass eigenstates 



V=2pJ2Pn [\(D + \H w \n)\ 2 - \(D-\H w \n)f 



(10) 



A similar pair of expressions can be written for x, 



X = f 



1 
f 



(D°\H W \D°)+PJ2 

n 

(d°\h w \d°)+py; 



(D°\n w \n){n\n w \D°) + (D°\H w \n)(n\H w \D°) 



m 2 D -El 



\(D + \H w \n)\ 2 -\(D_\H w \n)\< 
m 2 n - El 



(11) 



Note that x and y are generated by off-shell and on-shell intermediate states, respectively. 



III. SU(3) ANALYSIS OF D° - D° MIXING 

We now prove that D° — D° mixing arises only at second order in SU(3) breaking effects. 
The proof is valid when SU (3) violation enters perturbatively. This would not be the case, 
for example, if D transitions were dominated by intermediate states or single resonances 
close to threshold. As we will see explicitly in Sees. and 0, in such cases it is sometimes 



possible for SU(3) violation to be enhanced substantially. Yet other than in these exceptional 
situations, treating SU(3) violation perturbatively seems to us to be a mild assumption. 

The quantities Myi and Ti2 which determine x and y depend on matrix elements with 
the general structure 

74,74, ID ), (12) 
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where in this section we let 7i w denote specifically the AC = — 1 part of the weak Hamilto- 
nian. Let D be the field operator that creates a D° meson and annihilates a D°. Then the 
matrix element may be written as 

(0\DH w H w D\0). (13) 

Let us focus on the SU (3) flavor group theory properties of this expression. 

Since the operator D is of the form cu, it transforms in the fundamental representation 
of SU (3), which we will represent with a lower index, D^. We use a convention in which the 
correspondence between matrix indices and quark flavors is (1,2,3) = (u,d,s). The only 
nonzero element of Di is D\ = 1. The AC = — 1 part of the weak Hamiltonian has the flavor 
structure (qic)(qjqk), so its matrix representation is written with a fundamental index and 
two antifundamentals, H l k 3 . This operator is a sum of irreducible representations contained 
in the product 3x3x3 = 15 + 6 + 3 + 3. In the limit in which the third generation is 
neglected, H k is traceless, so only the 15 (symmetric on i and j) and 6 (antisymmetric on 
i and j) representations appear. That is, the AC = — 1 part of Ti w may be decomposed as 
±(015 + 6 ), where 

0jTj = (sc)(ud) + (uc)(sd) + si(dc)(ud) + Si(uc)(dd) 

— si(sc)(us) — si(uc)(ss) — s1(dc)(us) — sl(uc)(ds) , 
6 = (sc)(ud) — (uc)(sd) + si(dc)(ud) — si(uc)(dd) 

— si(sc)(us) + Si(uc)(ss) — Si(dc)(us) + sl(uc)(ds) , (14) 

and si = sm9 c ~ 0.22. The matrix representations H (15)^ and H(6)^ have nonzero 
elements 

If (IS)?' : = Hf = 1 , Hf = Hf = 8l 



H? = Hf = - 8l , HP = Hf = -4 , 

tt13 _ _ rj31 _ i rrl2 _ _ rr21 



(15) 



(6)Sf : H? = -Hf = 1 , Hf = -Hf = 8l , 

Hf = -Hf = - 8l , Hf = -Hf = -si . 

We introduce SU(3) breaking through the quark mass operator At, whose matrix represen- 
tation is Mj = diag(m u , ma, m s ). Although Ai is a linear combination of the adjoint and 
singlet representations, only the 8 induces SU (3) violating effects. It is convenient to set 
m u = m d = and let m s ^ be the only SU (3) violating parameter. All nonzero matrix 
elements built out of D i: H l k J and Mj must be 577(3) singlets. 

We now prove that D° — D° mixing arises only at second order in SU (3) violation, by 
which we mean second order in m s . First, we note that the pair of D operators is symmetric, 
and so the product DiDj transforms as a 6 under S77(3). Second, the pair of 7Y w 's is also 
symmetric, and the product H l k 3 H l ™ is in one of the representations which appears in the 
product 

"(15 + 6) x (15 + 6)1 = (T5 x T5) 5 + (15 x 6) + (6 x 6) s (16) 

-I s 

= (60 + 24 + 15 + 15' + 6) + (42 + 24 + 15 + 6 + 3) + (15' + 6) . 

A straightforward computation shows that only three of these representations actually ap- 
pear in the decomposition of H W H W . They are the 60, the 42, and the 15' (actually twice, 
but with the same nonzero elements both times). So we have product operators of the form 

DD = V 6 , 

'Hw'Hw = 060 + ^42 + 015', (17) 
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where the subscript denotes the representation of SU(3). 

Since there is no 6 in the decomposition of H W H W , there is no SU(3) singlet which can 
be made with and no SU(3) invariant matrix element of the form ([T3|) can be formed. 
This is the well known result that D° — D° mixing is prohibited by SU(3) symmetry. 

Now consider a single insertion of the SU(3) violating spurion M.. The combination 
T>§M. transforms as 6 x 8 = 24 + 15 + 6 + 3. Note that there is still no invariant to be made 
with H W H W . It follows that D° -D° mixing is not induced at first order in SU (3) breaking. 

With two insertions of Ai, it becomes possible to make an 577(3) invariant. The decom- 
position of VA4M. is 

6 x (8 x 8) 5 = 6 x (27 + 8 + 1) 

= (60 + 42 + 24 + 15 + 15' + 6) + (24 + 15 + 6 + 3) + 6. (18) 

There are three elements of the 6 x 27 part which can give invariants with 7i w 7i w . Each 
invariant yields a contribution proportional to s\m 2 . As promised, D° — D° mixing arises 
only at second order in the SU(3) violating parameter m s . 



IV. ESTIMATING THE SIZE OF SU(3) BREAKING 

We now turn to review some general estimates of the size of SU(3) breaking effects. 
These effects can be approached from either an inclusive or an exclusive point of view. It is 
instructive to see how SU (3) violation appears in each case. 



A. "Inclusive" approach 

An elegant and concrete estimate of how SU(3) violation enters x and y is the short 
distance analysis, first applied to D° — D° mixing by Georgi Q and later extended by other 
authors || |i~8f . We review it briefly, both to establish the contrast with our approach and 
to recall the results. Let A be a scale characteristic of the strong interactions, such as 
mp or Anf-x. In the limit m c ^> A, the momentum flowing through the light degrees of 
freedom in the intermediate state is large and an operator production expansion (OPE) can 
be performed. For example, one can write 

T 12 = ^- Im (D°\ i J d'xT{Ht C= \x) K C= \V)}\D") , (19) 

where H,^ c=1 is the |AC| = 1 effective Hamiltonian. In the OPE, the time ordered product 
in Eq. ([19]) can be expanded in local operators of increasing dimension; the higher dimension 
operators are suppressed by powers of A/m c . 

The leading contribution comes from the dimension-6 |AC| = 2 four-quark operators 
corresponding to the short distance box diagram, 

02 = u a 'j fl P L c l3 u l3 'ypP L c a , 0' 2 = u a P L cpUpP L c ai (20) 

where Pl — |(1 — 75). If one neglects QCD running between My/ and m c , in which case O2 
and 0' 2 do not contribute, one finds the simple expressions 



* n r 2 ^ (ml - mf,) 2 

AM box = — X D y 



3tt 2 " m\ 



5B' D m 2 D 



4 B D (m c + m u f 



(21) 
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ratio 


4-quark 6-quark 8-quark 


AM/AM box 


1 A 2 /m s m c {a s /Air)(A 2 /m s m c ) 2 


AT /AM 


m 2 /m 2 a s /&ir (3oa s /4ir 



TABLE I: The enhancement of AM and Ar relative to the box diagram at various orders in the 
OPE. A denotes a hadronic scale around 4^/^ ~ 1 GeV. 



Ar 



box 



4 (m 2 — m 2 \ 2 

— x D ^- 



miy m 2 s + m\ 



37r " m 2 m 2 



1 _ I?* — 

2 B D (m c 



7U 



D 



(22) 



where X D = V 2 s V 2 d G 2 F m D B D fp, and Bp are bag factors for , normalized to one in 
vacuum saturation. Including leading logarithmic QCD effects enhances this estimate of 
Ar by approximately a factor of two . Eqs. fl2"T|) and (^) then lead to the estimates 



•^bo 



few x 10 



Vbo 



few x 10" 



(23) 



Neglecting md/m s , Eq. (|22|) is proportional to m®. This factor comes from three sources: 
(i) m 2 from an SU(3) violating mass insertion on each quark line in the box graph 



ii 



m 2 from an additional mass insertion on each line to compensate the chirality flip from 
the first insertion; (iii) m 2 to lift the helicity suppression for the decay of a scalar meson 



into a massless fermion pair. The last factor of m 2 is absent from Eq. fl2~T|) for AM; this 
is why at leading order in the OPE ?/box "C a;b ox - Higher order terms in the OPE are 
important, because the chiral suppressions can be lifted by quark condensates instead of 
by mass insertions, allowing AM and Ar to be proportional to m 2 . This is the minimal 
suppression required by SU(3) symmetry, as we proved in Sec. |IT. 

The order of magnitudes of the resulting contributions are summarized in Table |. In 
the first line, the contributions to AM are normalized to AMb ox ; m the second line, the 
contributions to AT are normalized to AM at each order. The contribution of 6-quark 
operators to AM is enhanced compared to the 4-quark operators by A 2 /m c m s . This can 



be as much as an order of magnitude, if we identify the hadronic scale A as 47r/ T [25 



The second chiral suppression can also be lifted, but only at the price of adding a hard 
gluon, so the contribution of 8-quark operators to AM compared to the 6-quark operators is 
(a s /47i)(A 2 /m c m s ), which is of order unity. 1 In the case of AT, higher dimension operators 
are even more important [18 . A 6-quark operator, including a hard gluon to give an on-shell 



intermediate state, lifts both a chiral suppression and the helicity suppression. The 8-quark 
operators require a second intermediate particle to contribute to AT, which can be obtained 
by splitting the gluon already present for AM into a quark pair [I8[ , only costing a factor of 



/3oa s /(4:iT) ~ 1, where (3q — 11 — -^rif = 9. Thus, the dominant contributions to x are from 
6- and 8-quark operators, while the dominant contribution to y is from 8-quark operators. 
With some assumptions about the hadronic matrix elements, the resulting estimates are 



x 



y 



< 10" 



(24) 



1 We disagree with Ref. [Q, in which it was claimed that x and y can arise at first order in m s . Such 
contributions were claimed to come from pseudogoldstone loops which diverge in the infrared. However, 
there are no such divergences because the n, K and r\ are coupled derivatively. Such a contribution would 
also be in conflict with our proof in Sec. [II that D mixing is second order in SU(3) violating effects. 
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It is a general feature of OPE based analyses that x ^ y. We emphasize that at this time 
these methods are useful for understanding the order of magnitude of x and y, but not for 
obtaining reliable quantitative results. For example, to turn the estimates presented here 
into a systematic computation of x and y would require the calculation of almost two dozen 
nonperturbative matrix elements. 



B. "Exclusive" approach 

A long distance analysis of D mixing is complementary to the OPE. Instead of assuming 
that the D meson is heavy enough for duality to hold between the partonic rate and the sum 
over hadronic final states, here one assumes that D transitions are dominated by a small 
number of exclusive processes, which are examined explicitly. This is particularly interesting 
for studying Ar, which depends on real final states in D decays. 

For a long distance analysis, it is useful to express the width difference directly in terms 
of observable decay rates. From Eq. (0), we find 



V = E^ckmH Vcp(n) cos5 n y B(D° -> n) B(D° -> n) , (25) 

n 

where 8 n is the strong phase difference between the D° — > n and D° — > n amplitudes. In 
decays to many-body final states, the strong phases may have different values in different 
regions of the Dalitz plot, in which case the sum is supplemented by an integral over the 
Dalitz plot for each final state. The CKM factor is ?7ckm — ( — l) ns > where n s is the number of 
s and s quarks in the final state. For example, r] CKM (K + K~) = +1 and r] CKM (K + 'K^) — — 1. 
The factor rjcp = ±1 is determined by the CP transformation of the final state, CP\f) = 
Vcp\f)i which is well-defined since |/) and |/) are in the same SU(3) multiplet. This factor 
is the same for the whole multiplet. For example, t]cp = +1 for the decays to K + K~, 
and therefore to all decays into two pseudoscalars. For states where different partial waves 
contribute with different CP parities, rjcp is determined separately for each partial wave. 
For example, Vcp{p + p) — +1 for p + p~ in a relative s or d wave, and —1 in a p wave. 
Finally, it is convenient to assemble the final states into SU (3) multiplets and write 



V = ^2va, y a = Vcp(a) Vckm(ti) cos 5 n \J B(D° -> n) B(D° -> n) , (26) 



where a indexes complete SU(3) multiplets. By multiplets we refer to the SU(3) represen- 
tation of the entire final state, not of the individual mesons and baryons. 

In practice, we cannot use Eq. fl2~6|) to get a reliable estimate of y, since the doubly Cabibbo 
suppressed rates have large errors, and there are very little data on strong phase differences. 
To proceed further, we would be forced to introduce model dependent assumptions about 
the amplitudes and/or their strong phases. For example, in two-body D decays to charged 
pseudoscalars (7r + 7r~, ti + K~, K + n~, K + K~), the SU(3) violation can enter through the 
decay rates or the strong phase difference. We know experimentally that in some of these 
rates the SU(3) breaking is sizable; for example B(D° -> K+K~) /B(D° -> tt+tt") ~ 2.8 |6| . 



Such effects were the basis for the claim in Ref. |2T| that SU(3) is simply inapplicable to 
D decays. In contrast, we know very little about the strong phase 5 which vanishes in the 
SU(3) limit; Ref. [0 presented a model calculation resulting in cos 5 ^ 0.8, but it is also 



possible to obtain much larger values for 5 f2~2 |. Using Eq. ( p6|) , the value of y a corresponding 
to the [/-spin doublet of charged n and K is 



y vK = B(D° -> 7T+7T-) + £(77° -> 7T) - 2 cos 5 Jb(D° 7<~-7r+) £(77° 7i~+tt-) . 

(27) 

The experimental central values, allowing for D mixing in the doubly Cabibbo suppressed 
rates, yield y n K — (5.76 — 5.29 cos 5) x 10~ 3 f|. For small 5 there is an almost perfect 
cancellation even though the ratios of the individual rates significantly violate 577(3). In 
the "exclusive" approach, x is obtained from y by use of a dispersion relation, and one 
generally finds x ~ y. 

At this stage, one cannot use the exclusive approach to predict either x or y. Any 
estimate of their sizes depends on computing SU(3) breaking effects. While this problem 
is not tractable in general, one source of SU (3) breaking in y, from final state phase space, 
can be calculated with only minimal and reasonable assumptions. We will estimate these 
effects in the next section. 



V. SU(3) BREAKING FROM PHASE SPACE 



We now turn to the contributions to y from on-shell final states. There is a contribution 
to the D° width difference from every common decay product of D° and D°. In the 577(3) 
limit, these contributions cancel when one sums over complete SU(3) multiplets in the final 
state. The cancellations depend on SU(3) symmetry both in the decay matrix elements and 
in the final state phase space. While there are certainly SU(3) violating corrections to both 
of these, it is extremely difficult to compute the 577(3) violation in the matrix elements in a 
model independent manner. 2 However, with some mild assumptions about the momentum 
dependence of the matrix elements, the 577 (3) violation in the phase space depends only on 
the final particle masses and can be computed. In this section we estimate the contributions 
to y solely from SU(3) violation in the phase space. 3 We will find that this source of 577(3) 
violation can generate y of the order of a percent. 

The mixing parameter y may be written in terms of the matrix elements for common 
final states for D° and D° decays, 

|/=^E/ [ p - S -]« n w \n)(n\ H w |77°> , (28) 

where the sum is over distinct final states n and the integral is over the phase space for state 
n. Let us now perform the phase space integrals and restrict the sum to final states F which 
transform within a single SU(3) multiplet R. The result is a contribution to y of the form 




The SU(S) breaking in matrix elements may be modest even in cases such as D — > K + K and D — > tt + tt 



for which the ratio of measured rates appears to be very far from the SU(3) limit |28 . 
3 The phase space difference alone can explain the large SU (3) breaking between the measured D — » K*tv 
and D — > piv rates, assuming no SU(3) breaking in the form factors p9[ . Recently it was shown that the 
lifetime ratio of the D s and D° mesons may also be explained this way ]30|. 
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where p n is the phase space available to the state n. In the SU(3) limit, all the p n are the 
same for n e Fr, and the quantity in braces above is an SU(3) singlet. Since the p n depend 
only on the known masses of the particles in the state n, incorporating the true values of p n 
in the sum is a calculable source of 577(3) breaking. 

This method does not lead directly to a calculable contribution to y, because the matrix 
elements (n\H. w \D°) and (D°\7i w \n) are not known. However, CP symmetry, which in the 
Standard Model and almost all scenarios of new physics is to an excellent approximation 
conserved in D decays, relates (D°\7i w \n) to (D°\7i w \n) . Since |n) and \n) are in a common 
577(3) multiplet, they are determined by a single effective Hamiltonian. Hence the ratio 



Vf,r = 



ZneF R (D°\ 




n)p n (n\ 




D°) _ EneF R (D°\ 




n)p n (n\ 


7~(-w 


D°) 




7~tw 


n)p n (n\H w 


D°) EneF R T(D° - n) 



(30) 



is calculable, and represents the value which y would take if elements of Fr were the only 
channel open for D° decay. To get a true contribution to y, one must scale yF,R to the total 
branching ratio to all the states in F R . This is not trivial, since a given physical final state 
typically decomposes into a sum over more than one multiplet Fr. The numerator of yF,R 
is of order sf while the denominator is of order 1, so with large 577 (3) breaking in the phase 
space the natural size of yF,R is 5%. 

In this analysis, phase space is the only source of SU (3) violation which we will include. 
Of course, there are other SU (3) violating effects, such as in matrix elements and final state 
interaction phases. The purpose of our calculation is to explore the rough size of 577(3) 
violation in exclusive contributions to y. We assume that there is no cancellation with other 
sources of 577(3) breaking, or between the various multiplets which occur in D decay, that 
would reduce our result for y by an order of magnitude. This is equivalent to assuming that 
the D meson is not heavy enough that duality can be expected to enforce such cancellations. 

We begin by computing yp,R for D decays to states F = PP consisting of a pair of 
pseudoscalar mesons such as ir, K, 77. We neglect f] — i]' mixing throughout this analysis, and 
we have checked that this simplification has a negligible effect on the numerical results. Since 
PP is symmetric in the two mesons, it must transform as an element of (8 x 8)5 = 27 + 8 + 1. 
In principle, there are three possible amplitudes for D° — > PP, one with the pair in a 27 
and TC W in a 15, 

A 27 (PP 27 )^H^D m , (31) 
one with the pair in an 8 and and 7i w in a 15, 

AfiPP^HlW, , (32) 

and one with the pair in an 8 and and Tt w in a 6, 

AKPP^Dj. (33) 

However, the product H l k 3 Dj with (ij) symmetric (the 15) is proportional to H^Dj with 
(ij) antisymmetric (the 6), and the linear combination A$ = TLg 5 — A\ is the only one which 
appears. Thus there are effectively two invariant amplitudes. There is no SU(3) invariant 
amplitude to produce the final state in an singlet. Note that since we are assuming 577(3) 
symmetry in the matrix elements, such final states do not appear in our analysis. 
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It is straightforward to use these invariants in Eq. ([3(]) to compute Df,r- As an example, 
for yppfi we obtain 



ypp,8 



~ $(77, 7?) + ~ $(7T°, 7T°) + \ $(77, 7T°) + $(tT+, TT") + $(#+, IT") - J $ (r/, X°) 
2 2 3 6 

- 1 $(77, K°) - $(K + , TT") - $(#-, 7T+) - i $(i^°, tt ) - i $(Z°, vr°) 
6 2 2 

i $(?7, +*(K~, 7T+) + i $(F°, TT ) + 0{S\) 

6 z 



(34) 



where $(Pi, P2) is the phase space integral for the decay into mesons Pi and P 2 . In a two- 
body decay, $(Pi,P 2 ) is proportional to |p| 2 ^ +1 , where p and £ are the spatial momentum 
and orbital angular momentum of the final state particles. For D° — > PP, the decay is into 
an s wave. It is straightforward to compute the required ratios from the known pseudoscalar 
masses, 

ypP8 = -0.0038 si = -1.8 x 10~ 4 , y PP ,27 = -0.00071 s\ = -3.4 x 10~ 5 . (35) 

These effects are no larger than one finds in the inclusive analysis. This is not surprising, 
since as in the parton picture, the final states are far from threshold. 

Next we turn to final states of the form PV, consisting of a pseudoscalar and a vector 
meson. Note that three-body final states 3P can resonate through PV, and so are partially 
included here. In this case there is no symmetry between the mesons, so in principle all 
representations in the combination 8x8 = 27+10 + 10 + 85 + 8^ + 1 can appear. For 
simplicity, we take the quark content of the <fi and u respectively to be ss and (uu + dd) / a/2, 
and consider only the combination which appears in the SU(3) octet. We have checked that 
reasonable variations of the <fi — u mixing angle have a negligible effect on our numerical 
results. For each representation, there is a single invariant, up to the same degeneracy for 
the 8 as in the PP case. Along with the analogues of Eqs. (|3"1"|)-(|3"3"D with coefficients P27 
and B 8 = Pg 5 — Pf , we have the new invariants 

B w (PV 10 ) t3k H%D n e kmn (36) 

for H w in a 15, and 

B^PV^H^D^ (37) 

for 7i w in a 6. It turns out that these two invariants are proportional to each other. As 
before, the SU (3) singlet final state is not produced. 

Both because one of the particles is more massive, and because the decay is now into a 
p wave, the phase space dependence is stronger than for the PP final state. We obtain the 
ratios 

ypV8s = 0.031 si = 0.15 x 10~ 2 , y PV ,8 A = 0.032 sj = 0.15 x 10~ 2 , 
ypvm = 0.020 s{ = 0.10 x 10~ 2 , y PVj o = 0.016 s\ = 0.08 x 10~ 2 , 
ypV27 = 0.040 si = 0.19 x 10" 2 . (38) 

For any representation of the final state, the effects are less than one percent. 

For the VV final state, decays into s, p and d waves are all possible. Bose symmetry 
and the restriction to zero total angular momentum together imply that only the symmetric 
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Final state representation 
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PP 


8 
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—0.0038 
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27 


-0.00071 


-0.0034 


PV 


85 


0.031 


0.15 




8 A 


0.032 


0.15 




10 


0.020 


0.10 




TO 


r\ r\~\ a 


n no 
U.Uo 




27 


0.040 


0.19 


{VV) s - waV Q 


8 


n no 1 
— U.Uoi 


— U.oy 




27 


—0.061 


—0.30 


(VV)p-wave 


8 


-0.10 


-0.48 




27 


-0.14 


-0.70 


( VV )d-wave 


8 


0.51 


2.5 




27 


0.57 


2.8 



TABLE II: Values of Hf,r for two-body final states. This represents the value which y would take 
if elements of F R were the only channel open for D° decay. 



SU(3) combinations appear. Because some VV final states, such as 4>K*, lie near the D 
threshold, the inclusion of vector meson widths is quite important. Our model for the 
resonance line shape is a Lorentz invariant Breit-Wigner normalized on < m < 00, 

m 2 r 2 

/(m; m R , T R ) = N(m R , T R ) ^ _ -^JL _ , (39) 

where m R and T R are the mass and width of the vector meson, and m 2 is the square of its 
four-momentum in the decay. For s wave decays, we find the ratios 

yvV8 = -0.081 sj = -0.39 x 10" 2 , y vv ,27 = -0.061 s\ = -0.30 x 10~ 2 , (40) 

while for p wave decays we find 

yvvfi = -0.10 sj = -0.48 x 10" 2 , y VV:27 = -0.14 sj = -0.70 x 10" 2 , (41) 

and for d waves, 

yvV8 = 0.51 sj = 2.5 x 10~ 2 , y VVt 27 = 0.57 sj = 2.8 x 10" 2 . (42) 

With these heavier final states and with the higher partial waves, we see that effects at 
the level of a percent are quite generic. The vector meson widths turn out to be quite 
important; if they were neglected, the results in the p- and d-wave channels would be larger 
by approximately a factor of three. The finite widths soften the £77(3) breaking which 
otherwise would be induced by a sharp phase space boundary. We have checked that our 
results are not very sensitive to variations in the line shape used to model the vector meson 
widths. Again, AP and PPV final states can resonate through VV, so they are partially 
included here. Our results for two-body final states are summarized in Table ||. 

As we go to final states with more particles, the combinatoric possibilities begin to pro- 
liferate. We will consider the final states 3P and 4P, and for concreteness require that 
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the pseudoscalars be found in a totally symmetric 8 or 27 representation of SU(3). This 
assumption is convenient, because the phase space integration is much simpler if it can be 
performed symmetrically. These final states should be representative; we have no reason to 
believe that this choice selects final state multiplets for which phase space effects are par- 
ticularly enhanced or suppressed. Note that 3 x (15 + 6) contains no representation larger 
than a 27. 

In contrast to the two-body case, for three-body final states the momentum dependence 
of the matrix elements is no longer fixed by the conservation of angular momentum. The 
simplest assumption is to take a momentum independent matrix element, with all three final 
state particles in an s wave. In that case, we find 

y 3P8 = -0.48 sj = -2.3 x 1(T 2 , y 3Pi27 = -0.11 s? = -0.54 x 1(T 2 . (43) 

Note that the SU(3) violation is smaller for the larger multiplets, as more final states enter 
the sum. It may be that the 8 is in some sense an unusually small representation for three or 
more particles, and that this mode enhances the SU(3) violation by providing fewer distinct 
final states among which cancellations can occur. The enhancement of 2/3p,8 over y 3 p^7 is not 
a peculiarity of s wave decays. We have also considered other matrix elements; for example, 
if one of the mesons has angular momentum I = 1 in the D° rest frame (balanced by the 
combination of the other two mesons), then the ratios become 

y 3P8 = -1.13s? = -5.5 x 10~ 2 , y 3 p,27 = -0.074s 2 = -0.36 x 1(T 2 . (44) 

Alternatively, we could imagine introducing a mild "form factor suppression," with a weight 
such as n i7 y(l — m? /Q 2 ) _1 , where m 2 - = {pi + Pj) 2 , and Q = 2 GeV is a typical resonance 
mass. The result then changes to 

y 3P8 = -0.44 s\ = -2.1 x 1(T 2 , y 3P ,27 = -0.13 s 2 = -0.64 x 1(T 2 . (45) 

Finally, we have studied the final state with four pseudoscalars, with the mesons in an 
overall symmetric 8 or a symmetric 27. We take a momentum independent matrix element. 
There are actually two symmetric 27 representations; we call the 27 the representation of 
the form = [M^M™ M^M™ + symmetric — traces] and the 27' the one of the form 
flg = [M^M^MlMl + symmetric - traces]. Then we find 

2/4P, 8 = 3.3 s\ = 16 x 1(T 2 , y 4Pt27 = 2.2 s 2 = 11 x 1(T 2 . 

y 4P2r = 1.9 s \ = 9.2 x 1(T 2 . (46) 

Here the partial contributions to y are very large, of the order of 10%. This is not surprising, 
since 4P final states containing more than one strange particle are close to D threshold, and 
the ones containing no pions are kinematically inaccessible. There is no reason to expect 
SU(3) cancellations to persist effectively in this regime. Our results for 3P and 4P final 
states are summarized in Table fTJ. 

Formally, one could construct y from the individual yF,R by weighting them by their D° 
branching ratios, 

y=^E VfJ E r(D°-n)l. (47) 

L n£F R 1 



r 

1 F,R 



However, the data on D decays are neither abundant nor precise enough to disentangle the 
decays to the various 577(3) multiplets, especially for the three- and four-body final states. 
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Final state representation 
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-0.36 
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—0.44 
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27 


-0.13 


-0.64 


4P 


8 


3.3 


16 




27 


2.2 
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TABLE III: Values of vfr for three- and four-body final states. 



Final state 


fraction 


PP 


5% 


PV 


10% 


{VV) S - WSi yQ 


5% 




5% 


3P 


5% 


4P 


10% 



TABLE IV: Total D° branching fractions to classes of final states, rounded to nearest 



HI 



Nor have we computed Df,r for all or even most of the available representations. Instead, 
we can only estimate individual contributions to y by assuming that the representations 
for which we know Uf,r to be typical for final states with a given multiplicity, and then to 
scale to the total branching ratio to those final states. The total branching ratios of D° 
to two-, three- and four-body final states can be extracted from Ref. [p6[] . The results are 
presented in Table [TV], where we round to the nearest 5% to emphasize the uncertainties in 
these numbers. Close to half of all D° decays are accounted for in this table; the rest are 
decays to other modes such as PPV, decays to states with SU(3) singlet mesons, decays 
to higher resonances, semileptonic decays, and other suppressed processes. Based on data 
in the channel K°*p°, the VV final state is dominantly CP even, consistent with an equal 
distribution between s and d wave decays (although favoring a small s wave enhancement). 

We estimate the contribution to y from a given type of final state by taking the product 
of the typical yF,R found in our calculation with the approximate branching ratios given in 
Table |IV]. Such estimates are necessarily crude, but they are sufficient to give a sense of 
the order of magnitude of y which is to be expected. While in most cases the contributions 
are small, of the order of 10~ 3 or less, we observe that D° decays to nonresonant 4P states 
naturally contribute to y at the percent level. The reason for such unusually large SU (3) 
violating effects in y is that approximately 10% of D° decays are to final states for which 
the complete SU (3) multiplets are not kinematically accessible. 

It should be noted that for D decays to final states so close to threshold, our argument 
that D mixing is second order in SU(3) violation is inapplicable, because its underlying 
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assumption that SU(3) violation enters perturbatively is not met. In particular the proof 
fails near D threshold, if the decay is either to weakly decaying final states or to hadrons 
with widths T which are smaller than m s . In either case, the phase space available for 
the decay can vary rapidly on the scale of m s , spoiling the analytic expansion. For decays 
to hadronic resonances, T/m s is a small parameter which is not analytic as m s — > 0. For 
decays to long lived mesons, the O-functions which fix the boundaries of phase space are 
not analytic functions of their arguments, which in turn depend on m s through the masses 
of the final state hadrons. In this way, the generic m 2 s /m 2 D suppression is lifted and we find 
larger SU(3) violation in y just at the point that the conditions of the proof are not satisfied. 
We will see a similar failure of SU(3) cancellations when we study D mixing induced by 



resonances in Sec. VI 



We have not considered all possible final states which might give large contributions to 
y. In particular, the branching ratio for D° — > K~af is (7.3 ± 1.1)% [p6fl , even though this 
final state is quite close to D threshold. Unfortunately, the identities of the SU (3) partners 
of the ax (1260), which has J PC = 1++, are not well established. While it is natural to 
identify the ld(1400) as the corresponding strange axial vector meson, and the /i(1285) 
as the analogue of the u, there is no natural candidate for the ss analogue of the 0. The 
size of ypv* is quite sensitive to this choice, as well as to the value taken for the poorly 
measured width of the ai. If we take the ss state to be the /i(1420), and r(ai) = 400 MeV, 
we find ypv*,8 s = 1-8%. If instead we take the /i(1510), we find ypv*,8 s = 1-7%. With 
r(di) = 250 MeV, these numbers become 2.5% and 2.4%, respectively. Although it is clear 
that percent level contributions to y are possible from SU (3) violation in this channel, the 
data are still too poor to draw firm conclusions. 

On the basis of this analysis, in particular as applied to the 4P final state, we would 
conclude that y on the order of a percent would be completely natural. Anything an order 
of magnitude smaller would require significant cancellations which do not appear naturally 
in this framework. Cancellations would be expected only if they were enforced by the OPE, 
that is, if the charm quark were heavy enough that the "inclusive" approach were applicable. 
The hypothesis underlying the present analysis is that this is not the case. 



VI. SU(3) BREAKING FROM NEARBY RESONANCES 

One interesting feature of the D° is that there are excited mesons with masses close to 
mp>. As a result, it would not be unnatural for K resonances to play an important role in 
D decays. This possibility has already been explored in the literature |^U|, [27], |3T], ^] . Here 



we explore SU(3) breaking in the resonance contribution to D mixing. 

We are interested in the process D° — > R — > D°, where R is a resonance with mass mp 
and width Tp. Only spin zero resonances are relevant. The contribution of a single state to 
the D mass and width differences is given by 

res = \ H R? res = 2 \ H R? m l ~ ™R 

yR lR T {m 2 D - m 2 R f + m 2 D T 2 R ' R lR T m D (m 2 D - m 2 R Y + m 2 D T 2 R { ' 

where \H R \ 2 = {D°\Hw\R) {R\Hw\D Q ) parameterizes the couplings of R to D° and D°, and 
t]r is the CP eigenvalue of the 577(3) multiplet within which the resonance resides. If we 
assume the absence of direct CP violation in D decays, then (D \Hw\R) ma Y be related to 
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(R\Hw\D°) by SU(3). The ratio 

2j? _2{m 2 D - m 2 R ) 



y V K S m D T 



(49) 



is independent of Hr. Significant contributions to x and y from the resonance R are possible 
only if m 2 D — m\ < m D Y R . 

As a concrete example, consider the i\~*(1950), a positive parity excited kaon which, 
because of its large width, may play an important role in mediating D° — > K~tt + . Fitting the 
if* (1950) contribution to the observed D — > rates, one finds that resonance mediation 
could be as large as the usual quark tree amplitude P7|| . We can estimate an upper bound on 
the contribution of if* (1950) to y by assuming that the resonance is completely responsible 
for D —>■ Kir. The limit is given by 



(D° 


Hw 


K H )(K H 


Hw 


D°) 


(D° 


Hw 


K H ){K H 


Hw 


D°) 



^ /nnioy VH? lay mm X 15717 ~~r7Z\ » ( 5 °) 



I ATI 



where we denote the AT*(1950) by K H . With B(D° -> Kn) ~ 6% and i3(^ -> /Ctt) ~ 52%, 
we find |y| < 0.06 sj~3x 10 -3 . If D mixing is mediated by a resonance, then we expect x 
and y to be roughly of the same size. 

This upper bound is too generous, because we have not included the suppression from 



SU (3) cancellations. Note that our proof of Sec. |T|, that SU (3) violation appears only at 
second order in m s , applies only so long as m s Tr. While this must be true in the limit 
tur ~ mo —> oo, in which case Tr scales as m c , the ratio m s /r^ may not be small for 
resonances near the physical D mass. Therefore, SU(3) cancellations may be less effective 
for resonances than for real final states. 

The resonances in question fall into a positive parity 8 of 577(3), consisting of states 
which we will denote (n H , K H ,r] H ). If these states were degenerate and had equal widths, 
their contributions to D mixing would cancel. A measure of the actual effectiveness of this 
cancellation is the contribution of the entire multiplet relative to that of the Kh- The 577(3) 
partners of the TT*(1950) have not been conclusively identified. Instead of speculating, we 
will explore the efficiency of S77(3) cancellations qualitatively by taking the simple mass 
relations ^ 

m 7TH = m K[I - m s , m m = m K[I + - m s , (51) 
and assuming that the widths of the tth and tjh are the same as T(Kh) — 200 MeV. Then 

V-- 

4 yw H 4 



For m s = 150 MeV, we find y res /y T K S H — 0.27. The cancellations are somewhat less effective 
for x ves , with x res /x r £ s H = 0.50. We see that even for the if* (1950), likely to be the most 
favorable for inducing a large effect, S77(3) cancellations reduce the contributions to x res 
and y ies . We conclude that it would be quite unlikely for resonances to make a contribution 
to y at the level of one percent. 



VII. CONCLUSIONS 

The motivation most often cited in searches for D°—D° mixing is the possibility of observ- 
ing a signal from new physics which may dominate over the Standard Model contribution. 
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But to look for new physics in this way, one must be confident that the Standard Model 
prediction does not already saturate the experimental bound. Previous analyses based on 
short distance expansions have consistently found x, y ^ 1CT 3 , while naive estimates based 
on known SU(3) breaking in charm decays allow an effect an order of magnitude larger. 
Since current experimental sensitivity is at the level of a few percent, the difference is quite 
important. 

In this paper we have performed a general SU(3) analysis of the contributions to y. We 
proved that if SU (3) violation may be treated perturbatively, then D° — D° mixing in the 
Standard Model is generated only at second order in SU(3) breaking effects. Within the 
exclusive approach, we identified an SU(3) breaking effect, SU(3) violation in final state 
phase space, which can be calculated with minimal model dependence. We found that phase 
space effects alone can provide enough SU(3) breaking to induce y ~ 10~ 2 . Large effects 
in y appear for decays to final states close to D threshold, where an analytic expansion in 
SU (3) violation is no longer possible. 

We believe that this is an important result. Despite the large uncertainties, this is the first 
model independent calculation to give y close to the present experimental bounds. While 
some degree of cancellation is possible between different multiplets, as would be expected in 
the m c — > oo limit, or between SU (3) breaking in phase space and in matrix elements, it is 
not known how effective these cancellations are. The most reasonable assumption in light 
of our analysis is that they are not significant enough to result in an order of magnitude 
suppression of y. Therefore, any future discovery of a D meson width difference should not 
by itself be interpreted as an indication of the breakdown of the Standard Model. 

However, our analysis does not amount to a Standard Model calculation of y. First, we 
have considered only SU{3) breaking from phase space, and have not included any symmetry 
breaking in the matrix elements. Second, we have not calculated the contributions from all 
final states. Had we done so, we would still need very precise experimental data in order to 
disentangle the various SU (3) multiplets and combine the results into an overall value of y. 
Third, we have assumed that the charm quark is not heavy enough for duality to enforce 
significant cancellations between the various nonleptonic D decay channels, although some 
degree of cancellation is to be expected. 

The implication of our results for the Standard Model prediction for x is less apparent. 
While analyses based on the "inclusive" approach generally yield x ^ y, it is not clear 
what the "exclusive" approach predicts. The effect of SU (3) breaking in phase space in x 
is softer than in y, so one would expect x < y from our analysis. Thus if x > y is found 
experimentally, it may still be an indication of a new physics contribution to x, even if y is 
also large. On the other hand, if y > x then it will be hard to find signals of new physics, even 
if such contributions dominate AM. The linear sensitivity to new physics in the analysis of 
the time dependence of D° — > K + n~ is from x' = xcosS + ysin5 and y' = ycosS — :rsin<5 
instead of x and y. If y > x, then 5 would have to be known precisely for these terms to be 
sensitive to new physics in x. 

There remain large uncertainties in the Standard Model predictions of x and y, and 
values near the current experimental bounds cannot be ruled out. Therefore, it will be 
difficult to find a clear indication of physics beyond the Standard Model in D° — D° mixing 
measurements. We believe that at this stage the only robust potential signal of new physics 
in D° — D° mixing is CP violation, for which the Standard Model prediction is very small. 
Unfortunately, if y is larger or much larger than x, then the observable CP violation in 
D°-D° mixing is necessarily small, even if new physics dominates x. Therefore, searching 
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for new physics and CP violation in D° — D° mixing should aim at precise measurements 
of both x and y, and at more complicated analyses involving the extraction of the strong 
phase in the time dependence of doubly Cabibbo suppressed decays. 
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